Coherent effects in double-barrier Josephson junctions 
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The general solution for ballistic electronic transport through double-barrier Josephson junctions is 
derived. We show the existence of a regime of phase-coherent transport in which the supercurrent 
is proportional to the single barrier transparency and the way in which this coherence is destroyed 
for increasing interlayer thickness. The quasiparticle dc current at arbitrary voltage is determined. 
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Phase-coherent electronic transport in mesoscopic 
structures between normal (N) and superconducting (S) 
metals received considerable interest both in experiments 
and in theorya. Particular interesting phenomena were 
discovered in structures containing tunnel barriers (I). It 
is well known that the subgap resistance of a ballistic 
SIN junction has a quadratic dependence on the trans- 
parency of the interfacen, since Andreev reflection is a 
two-particle process. Disorder in the normal region en- 
hances the Andreev current due to opening of some frac- 
tion of tunneling channels, and the resistance has a linear 
dependence on the transparency. This effect isJmown as 
renectionless tunneling in SIN junctions (see&Q and fur- 
ther references iru). Interestingly, the opening of tunnel 
channels may be realized in a ballistic NI1MI2S junction 
as well by placing a second tunnel barriera. 

At the same time, a supercurrent I c in a tunnel SIS 
contact depends linearly on the barrier transparency 
since Cooper pairs tunnel coherentlyOu. Here we address 
the problem of universal features of supercurrent flow in 
a double-barrier ballistic SIiS'^S junction, where 5" is a 
thin layer with critical temperature T cs i < T cs . Coherent 
effects in such structures are also of practical importance, 
since recent experiments demonstrated the possibility of 
engineering Josephson junctions withpdesired properties 
using existing multilayered techniquea3-EI. 

The supercurrent in a disordered double-barrier 
SINIS junction was calculated by Kupriyanov and 
Lukichevli9, who considered the interlayer in the dirty 
limit and T cs > = 0. The coherent regime was found in 
the limit of small interlayer thickness d, when supercur- 
rent is of the first order in the single barrier transmission 
coefficient D, like in a SIS junction. However in the limit 
of small d the assumption of the dirty limit is not justified 
since electronic scattering takes place only at the inter- 
faces. Theoretical work on ballistic SINIS structures 
was concentrated on studying—resonant supercurrents in 
low-dimensional structuresO~li3. 

In the present paper we study theoretically the univer- 
sal features of charge transport in a three-dimensional 
(3D) SIiS'^S junction in the clean limit. We demon- 
strate the existence of the coherent regime when the su- 
percurrent is averaged over the transmission resonances 
and is proportional to D, whereas it becomes incoher- 
ent, of the order of D 2 , with increasing thickness, as 



expected for two uncorrelated sequential tunneling pro- 
cesses. We study quantitatively the crossover between 
these two regimes and the relation to the dirty limit re- 
sults of Ref.tx Further, we show that the coherent su- 
percurrent can be exactly derived from the distribution 
of transmission eigenvalues p(D) cx Z?~ 3 / 2 ll — D)~ x l 2 
known for a two-barrier NI1NI2N contactEl. Based on 
this distribution, we calculate the quasiparticle dc cur- 
rent at arbitrary voltage, which shows signatures of mul- 
tiple Andreev reflections (MAR). 

We consider a 3D ballistic SIiS'^S contact, where S' 
is a thin superconducting film with T cs i < T cs and the 
mean free path l s i » d, d is the interlayer thickness 
and I\ y 2 are the parallel atomically sharp interfaces with 
arbitrary transmission coefficients. In the temperature 
Green's function method the supercurrent density J s is 
expressed through the Fourier transform of the Green's 
function G (r, r') over the transverse coordinates 
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where x, x' are the coordinates across the junction, k\\ is 
the wave-vector component in the junction plane, u) n = 
(2n + 1)ttT . The normal and the anomalous Green's 
functions G(x, x'), F + (x, x') obey the Gor'kov equations 
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where A = A exp(ix) is the spatially dependent com- 



plex pair potential, H = -i— 



2m dx 2 
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V{x), E x = 

Ef — Ti z k 2 /2m is the electron kinetic energy across the 
junction, Ep is the Fermi energy, V(x) = W\5 (x) + 
W28 (x — d) is the interface potential, W\ 2 being the bar- 
rier strengths. 

Let us choose the position x' within the interlayer, 
then the solution of Eq. (||) in the superconducting 
electrodes (\x\ > d/2) is given by a linear combina- 
tion ofplane waves A(x') exp(ik x x) . Substituting it into 
Eq. (|2|), we arrive the dispersion relation for k which 

yields four solutions: k = ± 
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2mE F /h 2 . As a result the homogeneous solutions of Eq. 



at x < —d/2 has the form 
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where /? = iA*/(u> n + £). All four terms in Eq. (0) 
decay properly at x = — oo. A similar solution holds for 

x > -d/2 
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In the case of T cs / = the solution in the interlayer 
\x\ < d/2 consists of the homogeneous part, in which f3 = 
and all terms of Eqs. ([| |^) are present, and the par- 
ticular solution (source term): G = {m/iTi 2 k)e lk \ x ~ x I, 
F = 0. 

The solutions in all three regions are matched by the 
conditions of continuity of the G(x,x') and F(x,x') at 
the interfaces x = ±d/2 and by the condition for the 
derivatives which follows from the integration of Eq. (j|) 
across the interface barriers. This yields for x = —d/2 
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and similar condition for x = d/2. These boundary condi- 
tions provide the required number of linear equations for 
the coefficients in the equations for G(x, x') and F(x, x'). 
We can expand the wavevector in the interlayer as 
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where k x = */ k F — fcjj = kp cos 6 is the transverse com- 
ponent of the wave vector, £ s , x — £ s , / cos6>, £ s , — 
Tivf/2 cos 0-^/w 2 + A 2 , is the coherence length and Ay 
is the pair potential in the S' layer. In the S electrodes 
we still keep k$ = k x with accuracy up to terms of the 
order of A/Ep. By matching the solutions |^) in all 
three regions at the interfaces we derived an expression 
for the supercurrent J s , valid for arbitrary d and W\^.. 
Below we present the solution for d < £ s 'and symmetric 
low-transparent barriers Wi^ = W, (W/hvp) ^> 1 
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A 2 sin ip + A s A fl . ^W/E~ 2 d/i sl W 2 sin | 
2 W 4 E 2 (cosh d/£ s , x - cos 2k x d) + E 2 
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Here p is the phase difference across the junction, A s 



A 2 , E 2 = 
2 + A 2 cos 2 ip/2 and W = W/hv F . 



is the pair potential in S, E\ 

Eq. (@) is the main technical result of this paper and 
describes the interplay between quasiparticle transmis- 
sion (Brcit-Wigner) resonances and electron-hole (An- 
dreev) resonances. Changing the phase space in inte- 
gration over fc|| , one can apply Eq.(||) to the prob- 
lem of the supercurrent flow via transmission resonances 
irjpJcpfcdimpnsional contacts. In this case the results 
oejcjO'Ej, taken in relevant limits, are reproduced. 
Below we concentrate on the interplay between coher- 
ent and incoherent regimes in the case of large k F d. 
It is instructive to discuss this crossover in terms of 
the width of transmission resonances in a double-barrier 
junction, which in the symmetrical case is given by 
T = hv F (xD(x)) )2d = hv F /8dW 2 , where (xD(x)) 
is the angle-averaged transparency of a single barrier 
(x = cos 9). 

Coherent regime (broad resonances). As follows from 
Eq.(|J), the coherent regime takes place for a thin inter- 
layer when the transmission resonances are broad T 3> 
7rT cs . The supercurrent is then given by 
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where D 2 is the transparency of a double barrier NININ 
contact 

D^ 1 = 1 + (ZW cos k'pd + 2W 2 sin k'pdj ^ , (10) 

and has a resonant structure. Integration over the direc- 
tions of fcii (over the resonances) yields the supercurrent 



J S R N = 2-kT 
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which does not depend on the properties of the intpt- 
layer and coincides with the dirty-limit KL resultlill 
Here Rjj 1 = e 2 kp^f/^-K 2 h is the normal state contact 

resistance per square, where 7 = (xD(x)} — 1/4W 2 . 
This expression can be generalized to the asymmet- 
ric case: eJ s RN — 27rT^ w l^s| 2 smp/EiE' 3 , with 



2 + |A s r(cos 2 ^/2 
72), 71,2 = (xDi, 2 (x)} 
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<P/ 2 ), 7- = (7i ~ 
where D\ 2 are the 



E' 3 = 
72)/ (71 

individual barrier transparencies. In this case, R^ 1 = 
e 2 kp'y c /2Tr 2 'h, where j c =^72/(71 + 72) which is the 
classical result known fromtj. 

Expression ( pr[ ) has been found by Kupriyanov and 
Lukichevlij in the case of a double-barrier junction with 
a dirty metal interlayer in the limit of vanishingly small 
thickness. This fact shows that Eq. ([ll]) is a very general 
result. 

For T — the maximum value of eI s RN is achieved at 
p as 1.86 and exceeds the eJ c RN value of (7r/2) |A s | for 
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FIG. 1. I c Rn product in the model of clean interlayer 
(solid lines) and disordered interlayer (dashes lines). The 
curves are plotted for the NB/A1 case for T cs — 7.4 T C3 i. 

a tunnel SIS contact. The reason is that in the coherent 
regime the dominant contribution to J s comes from the 
transmission resonances, which in the present case are 
broader than A. As a result the supercurrent is of the 
first order in (xD(x)). 

The supercurrent in the coherent regime has the spec- 
tral density 



Im J s (s) 
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for A s cos<^/2 < e < A s , while lmJs(e) = for e < 
A s cosip/2 and e > A s , i.e. the Andreev bound states in 
the energy range A s cos<p/2 < e < A s contribute to the 
supercurrent. 

Incoherent regime (narrow resonances) T <C 7rT cs . 
With the increase of the interlayer thickness the coherent 
regime breaks down due to the dephasing of the transmis- 
sion resonances. After performing the angle averaging in 
Eq.(||) the general expression for a double-barrier junc- 
tion becomes 
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2 i I (A 2 cos 2 f+^)/(A 2 +^). 
The pair potential in S' is determined self-consistently 
with 



A, In-^ = 27rTW^-<F s , 



(14) 



where (Fs>) is the angle-averaged anomalous Green's 
function, that is solved in the same way as the normal 
Green's function G. 
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As seen from Eq. (J13|) , the current-phase relation has 
two components, simp and sin<y9/2. The sign of the 
siny>/2 component is determined by the sign of the in- 
terlayer pair potential A s > , which is determined selfcon- 
sistently from Eq.(jlJ) and depends on the sign of the 
electron-electron interaction in an 5' material. For an 
attractive interaction (T cs > > 0) A s / does not vanish even 
at T > T cs ' and has a positive sign, while for a repulsive 
interaction A s / becomes negative. Therefore the mea- 
surements of a current-phase relation in a SIS' IS junc- 
tion can be used for measuring the sign of the electron- 
electron interaction in metallic films. 

As is shown above, the critical current is controlled 
by a single suppression parameter j e ff = ttT c /T — 
2TrT cs d/hv F (xD{x)). As follows from Eq.@, for T > 
T cs ' and j e ff(uj n /TtT cs ) 3> 1 the supercurrent becomes 
of the order of (xD(x)) 2 as expected for the incoherent 
tunneling in a double-barrier contact: 



eJ s RN — 



32-kT ^ A 2 sin ip 
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A numerical evaluation of Eq.(|i~3|) is shown in Fig. 1. 
The coherent (I) and incoherent (II) regimes are indi- 
cated. Regime number III shows the crossover to the 
series connection of two SIS' tunnel junctions. 

Transmission distribution, MAR. The result for the 
coherent regime Eq. (|lT|) can be derived from the trans- 
mission eig envalue density p{D) = (Gjv / 7 ipo)D~ 3 / 2 (l - 
D)^ 1 / 2 for two-barrier NININ contacts!!, where Go = 
e 2 /2ttTl. While deriving this distribution, an assumption 
about a certain amount of impurity scattering in the in- 
terlayer was madeEl; it can be shown that it also holds for 
the considered case of a clean interlayer, provided that 
k F d > 1. 

The derivation of Eq.(fTl"|) is then performed by the 
calculation of the integral I C (D) p(D)dD, where 
I C {D) is a supercurrent per single ballistic channel, 
IJD) = (e/2h)A 2 s Dsimpt&nh(E B /2T)/E B and E B = 

A\J \ — D sin 2 <p/2, by the residues of tanh(x) in a com- 
plex plane z = (1/ D— 1) x l 2 . It yields exactly the coherent 
supercurrent cJ s Rn — 27rT^ w>0 A 2 s\mp/ E\E%. This 
proof also holds for the asymmetric case, when the eigen- 
value density is p(D) = (G N /nG )D- 3 / 2 (D max - D)~^ 2 
with D max = 47172/(71 + I2) 2 , < D < L» max . 

The considerations above allow one to calculate the 
current under finite voltage V at the contact on 
the basis of the MAR (multiple Andreev reflections) 
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tailed theory applicable for arbitrary J e ff is the subject 
of further study. 

In conclusion, the general solution for the supercurrent 
in double-barrier Josephson junctions is presented and 
the crossover from the coherent to the incoherent regime 
for increasing interlayer thickness is studied in detail. 
The parameter-free calculation of quasiparticle current 
is done in the coherent regime. The results have implica- 
tions for transport in Josephson junctions and multilayers 
engineered with modern techniques. 

The authors thank M.Yu. Kupriyanov, K.K. Likharev, 
Yu.V. Nazarov, LP. Nevirkovets and H. Rogalla for many 
useful and stimulating discussions. 




FIG. 2. Dc current component in SIS'IS contact in the 
coherent regime, together with the single-channel curves for 
various transmissions D. 

formalism&0. Below we discuss the dc current com- 
ponent in a symmetric SIS'IS contact. The dc current 
per single channel with transmisison D is given byc3 

irhI q (D,V)/e = eVD- J de tanh JL(1 - |a | 2 ) (17) 



2Re(a A )+J2(\ A r l 



where a n — a(e + neV) is the Andreev reflection ampli- 
tude: a(e) = (s - sgn(e)(s 2 - A 2 ) 1 / 2 ) for |e| > A and 
a(e) = (e - i(A 2 - e 2 ) 1/2 ) for |e| < A. The coefficients 
A ni B n are-given by the recurrency relations which were 
derived inE3 and can easily be solved numerically. 

The resulting total dc current in a symmetric SIS' IS 
contact is given by the integration of the single- 
channel result I q (D, V) over the eigenvalue density I q — 
/o I q (D,V)p(D)dD where p(D) = (G N /ttG )^- 3/2 (1 - 
D)-Va. 

The results of numerical calculation at temperatures 
T <C T cs are presented in Fig. 2. A few single-channel 
curves for various D are also shown for comparison. 
The excess current I ex eRN — 1.05A S is present at high 
bias eV 3> A s , while the subharmonic gap structure at 
eV = 2A s /n due to MAR is present at lower voltages, 
despite averaging over the channels. This dc current de- 
termines the amount of dissipation in the junction. Note 
that the curve SIS'IS in Fig. 2 is universal, i.e. indepen- 
dent of microscopic parameters such as electronic mean 
free path or T c of the interlayer as long as j e ff < 1. This 
universality breaks down with the increase of Jeff, due 
to the dephasing of the transmission resonances. The de- 
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